Domain-wall structure of a classical Heisenberg ferromagnet on a Mobius strip 
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We study theoretically the structure of domain walls in ferromagnetic states on Mobius strips. 
A two-dimensional classical Heisenberg ferromagnet with single-site anisotropy is treated within 
a mean-field approximation by taking into account the boundary condition to realize the Mobius 
geometry. It is found that two types of domain walls can be formed, namely, parallel or perpendicular 
to the circumference, and that the relative stability of these domain walls is sensitive to the change 
in temperature and an applied magnetic field. The magnetization has a discontinuity as a function 
of temperature and the external field. 

PACS numbers: 75.60.Ch, 75.10.Hk, 75.60.Ej 



I. INTRODUCTION 

The effect of system geometry on the physical proper- 
ties has stimulated renowned interest after the reaUzation 
of crystals with unusual shapes, e.g., ring, cylinder and 
especially Mobius strip^ii^i^i^ because the ordered states 
in these systems could be quite different from those in the 
ordinary bulk systems. These single crystals are made of 
quasi-one-dimensional conductors, such as, NbSea. 

Hayashi and Ebisawa^ studied the s-wave supercon- 
ducting (SC) states on a Mobius strip based on the 
Ginzburg-Landau (GL) theory, and found that the Little- 
Parks oscillation is modified compared to that for the 
ordinary ring-shaped samples. When the flux quanta 
threading the Mobius strip are close to half-odd inte- 
ger times flux quantum (pg — hc/2e (h, c, and — e be- 
ing the Planck constant, the speed of light and the elec- 
tron charge, respectively), a state appears in which the 
real-space node in SC gap exists along the circumfer- 
ence of the strip. Later a more detailed examination us- 
ing the Bogoliubov-de Gennes (BdG) theory confirmed 
this prediction.^ Charge density wave (CDW) states in 
ring-shaped crystals were also examined by Nogawa and 
Nemoto,^ and Hayashi et al^ There the effect of frus- 
tration due to bending of the lattice was clarified. The 
persistent current in a Mobius strip as a function of the 
applied flux was examined by Yakubo et alJ^ and Mila 
et al.^^ Wakabayashi and Harigaya--^ studied the Mobius 
strip made of a nanographite ribbon, and the effect of 
Mobius geometry on the edge localized states has been 
explored. Kaneda and Okabe^ studied the Ising model 
on Mobius, especially the effect of sample geometry on 
the finite scaling properties. 

In this paper we study the ferromagnetic states on 
Mobius strips. Although the single crystals with Mobius 
geometry have been obtained for systems with CDW or 
superconducting order up to now, we expect that the 



synthesis of ferromagnetic Mobius strips using quasi-one- 
dimensional ferromagnets could be possible. In such a 
system there should be a domain wall (DW), and two 
kinds of DWs, i.e., those parallel and perpendicular to 
the circumference are expected. We will show that these 
two types of DWs can actually exist. An interesting point 
is that the relative stability of these two kinds of DWs is 
quite sensitive to the change in temperature and the ap- 
plied magnetic field. We expect that this might be used 
in the field of technological applications in future. 

This paper is organized as follows. In Sec. II the model 
and the method of calculations are presented. We study 
the structure of domain walls on Mobius strips in Sec. 
HI. In Sec. IV the effect of an applied magnetic field is 
examined. Section V is devoted to summary. 



II. MODEL AND MEAN-FIELD 
APPROXIMATION 

In order to study the ferromagnetic states on Mobius 
strips we consider a classical Heisenberg model with 
single-ion anisotropy, since this model is most suitable 
to examine the spatial variations of the magnetization. 
A model system with Mobius geometry is constructed 
as follows. We consider Ny (even integer) sets of ferro- 
magnetic chains (each chain has sites; being an 
integer), and these chains are weakly coupled ferromag- 
netically. We denote the i th spin on the j th chain 
(l < i < N^, 1 < j < Ny) as Sjj (|S,j| = 1). As a 
first step we form a cylinder, where the circumference is 
parallel to the chains. If one moves one site along the 
chain, the direction normal to the strip changes by the 
angle di = 2tt/Nx. Then we twist the strip along the cir- 
cumference at the middle of the strip (shown as a white 
curve in Fig.l) so that the twist angles at the neighboring 
sites on the white curve (for example, the points denoted 
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FIG. 1; Schematic picture of a Mobius strip. The axes are 
those for the relative coordinate system at each point. See 
text for details. 



as P and Q) differ by O2 = tt/N^- Finally the boundary 
condition is imposed, 



Sl,7V„ + l- 



(1) 



to get a Mobius strip. The direction normal to the strip is 
different from site to site, so we define the absolute coor- 
dinate system (x, y, z) and the relative coordinate system 
(i, y, z) at each point and use one of them which is better 
to describe the particular situation. For the former, the 
point with z = 1 on the middle of the strip {P in Fig.l) 
is taken to be the origin, and the z direction is perpen- 
dicular to the strip. The x and y directions are within 
the strip, and the former (latter) is parallel (perpendic- 
ular) to the circumference. For the relative coordinate 
system, the z direction is defined as that perpendicular 
to the strip at each point, and x (jj) direction is par- 
allel (perpendicular) to the circumference at that point 
(see Fig.l). Namely, the relative coordinate system for a 
particular point P is defined as the absolute coordinate 
system. Then the spin in the relative coordinate system 
at a site (j,j), S^j-, is expressed by that in the absolute 
coordinate system, S^j-, as 

Sf^ = cos(fc6li) - sm{k9i), 

st^ slcosike^) 

-{S^j sin(fc6li) + cos(fc6li)) sin(A:6i2), (2) 

St, - (Sf J cos(fc0i) + S't-sin{k6i)) cos(fc02) 

'v . 
-3 ■ 



sin(fc02), 



with k = i — 1. Now the Hamiltonian of our system is 
given as 



H — — J^^(Sij- • Si+i.j -I- aSij- • Sij+i) 



(3) 



We treat H. using a mean-field (MF) approximation 
and the resulting MF Hamiltonian reads 



(4) 



where va.i,j and Eq are the mean-field acting on the spin 
Si.j and the constant term of the energy, respectively. 



s=±i 



Eq^ J ^(Sjj) • ((Sj+5j) +a{S,j+s)), 

i,j S=±l 



(5) 



with the boundary condition (Si^o) = {Si,Ny+i) — 0. The 
spins Si J are expressed by the angles 0ij and (fiij ; S^j — 
sin 9ij cos , Sfj = sin 9ij sin (jj.j , S^^ = cos O^j , 



Z- - — Tr- -p~'^'^''J 
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(6) 



with (3 = 1/T. (T is the temperature and we use units 
kg = 1.) The free energy is calculated as 



F = £;o-r^iogz,,,. 



(7) 



i,3 



Equations (4)- (6) constitute the self-consistency equa- 
tions which we solve numerically by the method of it- 
eration. If we obtain several solutions to satisfy the self- 
consistency equations for the same set of parameters in 
the model, the solution with the lowest free energy is the 
true one. 

Because of the anisotropy in spin space due to the 
presence of the D term, Mermin- Wagner theoreroi^ does 
not apply to our system. Then the finite Tc would be 
expected and our mean-field treatment strictly in two- 
dimensional systems should be justified in a qualitative 
sense. In order to relate our model to real samples which 
have finite width perpendicular to x and y directions, we 
assume that there is a weak three dimensionality (a finite 
but thin width) and that the system is uniform along the 
third direction. These assumptions allow us to consider 
the realistic systems by use of the model presented above. 



where J(> 0) and Ja are the exchange energies between 
nearest-neighbor spins located along and perpendicular 
to the chains, respectively (0 < a < 1). D is the single- 
ion anisotropy energy which favors spin alignment per- 
pendicular to the strip if £> > and H is the external 
magnetic field. 



III. STRUCTURE OF DOMAIN WALL 



In this section we examine the spatial variations of the 
magnetization (S'j^j) (p = x,y,z) in order to clarify the 
structure of DWs on Mobius strips. In the following we 
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fix the system size to be = 100 and Ny = 10, and 
take J = 1 as the unite of energy. First we show the i 
and j dependences (direction parallel and perpendicular 
to the circumference, respectively) for D = 0.5, a = 0.04 
and T = 0.2Tc (Tc = 0.811J) in Fig.2. As seen in 
Figs. 2(a) and 2(b) (S^j) [fi = x, y, z; note that {x, y, z) is 
the absolute coordinate system] for i = in the j = 1 
chain connect smoothly to those for i — 1 in the j — 10 
chain. On the contrary in Fig. 2(c) there is a sign change 
in {S^ j) in the middle of the strip (between j = 5 and j = 
6) and this kind of sign change occurs for all i. It means 
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FIG. 3; Two types of domain walls: (a) type A - parallel 
to the circumference and (b) type B - perpendicular to it. 
{x,y, z) are the axes of the absolute coordinate system. 



pling, a = 0.15. In Fig. 4(a) the i dependence of {S^^) is 
shown for j = 1 and j = 10 for T = 0.2Tc (Tc = 0.883 J). 
(The results for j = 1 and j — 10 are the same.) In this 
case (S^j) for i = in the j = 1 chain do not con- 
nect smoothly to those for i = 1 in the j = 10, and 
so the discontinuity is observed. On the other hand the 
j dependence is smooth [Fig. 4(b)], and this behavior is 
the same for all i. It indicates that the DW is located 
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FIG. 2: The spatial variations of (Sij) for D = 0.5, a = 0.04 

and T = Q.2Tc. (a) j = 1, (b) j = 10 and (c) i = 51. Note FIG. 4: 

that {x, y, z) is the absolute coordinate system. and T - 



The spatial variations of {Si,j) for D 
-- 0.2Tc. (a) j = 1, and (b) i = 51. 



: 0.5, a = 0.15 



that a domain wall is formed along the circumference of 
the Mobius strip as schematically shown in Fig. 3(a), and 
we call this kind of DW as type-A . When we transform 
S'^j (absolute coordinate) to S^j (relative coordinate), 
the latter has only z component except near the DW. 
Next we consider the case of stronger inter-chain cou- 



perpendicularly to the circumference in contrast to the 
case of a = 0.04, and we call this kind of DW as type 
B [schematically shown in Fig. 3(b)]. Naively we expect 
that a type-A (B) domain wall is stable if N^a < Ny 
{Nr^a > Ny), since the energy to create a DW is smaller 
compared to that for type B (A). The results shown in 
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Figs. 2 and 4 are consistent with this expectation. 

However, this stability condition for DWs will not hold 
at higher temperatures or reduced values of D. In Fig. 5 
the phase diagram in the plane of a and T is shown. 
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FIG. 5: Phase diagram in the plane of T and a. Curves are 
the boundaries between type-A and B domain walls, and the 
region close to the origin corresponds to the type-A domain 
wall. 



As seen the region of type-B domain wall becomes larger 
as T is increased or D is decreased. This behavior can 
be understood as follows. Since we consider the quasi- 
one-dimensional system, the inter-chain exchange inter- 
action Ja is much smaller than the intra-chain exchange 
J. The single-ion anisotropy energy D was taken to be 
Ja < D < J, which seems to be a natural choice. Since 
the D term favors the spin alignment normal to the strip, 
the loss of energy would become large for the broader 
width of the DW of type A. Then the width of the type- 
A DW is quite small irrespective of the temperature, ac- 
tually a single site in most cases. On the contrary for 
the type-B DW, the width of the DW may become large 
because D < J. As D is further decreased or the temper- 
ature is increased, the width of DW can be larger and so 
the DW formation energy becomes smaller The spatial 
variations of {S"j ) (in the relative coordinate system) for 
85 < I < 100 and 1 < « < 15 are depicted in Fig.6. The 
width of the type-B DW is actually smaller at lower T 
for the same value of D. Here {S^j) for 85 < i < 100 
are multiplied by -1, since the normal direction in the 
relative coordinate system is reversed between the sites 
with i — 100 and i — I. We have also checked that the 
width becomes larger for smaller D at the same temper- 
ature, as stated above. Thus the type-B DW becomes 
favorable compared to type-A DW, as T is increased or 
D is decreased. 

From Fig. 5 we see that the phase transition from type- 
B DW to type-A DW occurs as the temperature T is 
lowered. This results in an interesting behavior of the 
magnetization as a function of T. In Fig. 7 the T depen- 
dences of the magnetizations for a ~ 0.15 and a = 0.04 
are shown, where the total magnetization for each di- 
rection, {fj, = x,y,z), and its magnitude, Mtot, are 
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FIG. 6: The spatial variations of (S^j) near the DW in the 
relative coordinate system for D — 0.5 and q = 0.15. T = 
O.lTc for (a) and (cj, and T = O.QTc for (b) and (d). 



defined as 



M„ 



" 1=1 J=l 



(8) 



In the case of a = 0.15 the magnetization monotonously 
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FIG. 7: T dependences of the magnetizations for D — 0.5, 
(a) a = 0.15, and (b) a — 0.04. The magnetization has a 
jump as indicated by the dashed vertical line. 



increases as T is lowered. On the contrary the magne- 
tization for a = 0.04 suddenly vanishes at T ~ 0.6J. 
This is because the phase transition from type-B DW 
to type-A DW occurs at this temperature. If the DW 
is of type A, the magnetizations in the j th chain are 
exactly canceled by those in the {Ny — j + 1) th chain, 
though locally there are finite magnetizations. On the 
other hand, in the case of type-B DW there is no such 
cancellation and the total magnetization is finite, and it 
becomes larger as T is lowered. In realistic systems this 
behavior of the magnetization would be observed with 
the hysteresis. However, the detection of this behavior 
may be difficult because type-A and type-B DW states 
have globally different structures so that the pinning of 
DWs due to impurities or dislocations would disturb their 
motion more seriously than usual cases. 

The Curie temperature for the system with Mobius 
geometry is only slightly lower than that for the planar 
system with the same parameters. For D = 0.5, a = 
0.04, = 100, and Ny = 10, the former is 0.811J while 
the latter is 0.812J. (The Tc for the infinite plane with 
the same D and a is 0.813J.) When and Ny become 
larger the difference of Tc will be smaller, since the local 
curvature of the Mobius strip is decreasing. This implies 
that the Curie temperature for the Mobius system is close 
to that for the flat system in general, at least in the mean- 
field approximation. 



IV. EFFECT OF EXTERNAL MAGNETIC 
FIELD 

Next we examine the effect of an external magnetic 
field H on the domain-wall structure. Here H is as- 
sumed to be applied along the y direction, i.e., H threads 
the Mobius strip. The T dependences of the magnetiza- 
tions are shown in Fig. 8. Because of finite H, the mag- 
netization is finite for T > Tc{H = 0) and gradually 
decreases as T is increased. For a ~ 0.15 [Fig. 8(a)] 




T 

FIG. 8: T dependences of the magnetizations for D — 0.5, 
H = 0.01, (a) a = 0.15, and (b) a = 0.04. The magnetization 
has a jump as indicated by the dashed vertical line. 

the magnetization smoothly grows as T is decreased for 
T < Tc{H = 0), while there is a discontinuity in the case 
of a = 0.04 [Fig. 8(b)]. This is again due to the transition 
from type-B to type-A DW, but the transition occurs at 
a much lower temperature compared to the case of if = 
[Fig. 7(b)]. This is because the Zeeman coupfing lowers 
the energy of the state with type-B DW due to its finite 
total magnetization and type-B state is more favored in 
the presence of the external field. For H the mag- 
netization in the type-A state is finite but is very small. 
Next we show the phase diagram in the plane of T and 
a for fixed value of H [Fig. 9 (a)] and that in the plane 
of T and H for fixed a [Fig.9(b)]. In both Figs.9(a) and 
9(b) the region of type-A DW shrinks when H is applied. 
This can be understood similarly as the T dependence of 
the magnetization. In the type-B state the magnetiza- 
tion is already finite without so that the energy can 
be gained due to the coupling to the applied magnetic 
field. Then the type-B state is more favorable compared 
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to the type-A state in the presence of H. This resulted 
in a serious difference in response to the appHed mag- 
netic field for two types of DWs. The dependence of the 
magnetization on the external magnetic field is shown in 
Fig. 10. There is again a jump of the magnetization as a 
function of H because of the transition between different 
types of DWs. If we assume J ~10meV, the transition 
point locates at around i?=lT for the parameters used 
here. 

In the above we have assumed that H is applied along 
the y direction. For the type-B DW, this direction is par- 
allel to the domain wall. When H is applied along other 
directions, type-B DW will move to a different position 
at which the more Zeeman energy can be gained. This 
problem will be examined elsewhere. 
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FIG. 9: Phase diagram (a) in the plane of a and T with D = 
0.5, and (b) in the plane of H and T with ol = 0.04. Curves are 
the boundaries between type-A and B domain walls, and the 
region close to the origin corresponds to the type-A domain 
wall. 
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In summary, we have studied the structure of domain 
walls in ferromagnetic Mobius strips. We found that 
there can be two kinds of domain walls whose relative sta- 
bility is quite sensitive to the change in temperature and 
the applied magnetic field. The magnetizations will have 
interesting behaviors in that they have discontinuities as 
functions of both temperature and the magnetic field. 
Although ferromagnetic Mobius strips have not yet been 
obtained, we expect that these systems will attract much 
interest both theoretical and experimental once they are 
synthesized. 
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